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We extend the original work of Ruderman, Kittel, Kasuya, and Yosida (RKKY) on the 
interaction between two magnetic moments embedded in an electron gas to the case where 
the electron gas is spin polarized. The broken symmetry of a host material introduces 
the Dzyaloshinsky-Moriya (DM) vector and tensor interaction terms, in addition to the 
standard RKKY term, so that the net interaction energy has the form: Pi = J S\ -S^ + D* 

Si x S 2 + Si- r -S 2 - We find that for the spin-polarized electron gas, a non-zero tensor 
^ —* 

interaction r is present in addition to the scalar RKKY interaction J, while D is zero 

due to the presence of inversion symmetry. Explicit expressions for these are derived for 
the electron gas both in 2D and 3D. The RKKY interaction exhibits a beating pattern, 
caused by the presence of the two Fermi momenta kp+ and kp\ , while the R~ 3 distance 
dependence of the original RKKY result for the 3D electron gas is retained. This model 
serves as a simple example of the magnetic interaction in systems with broken symmetry, 
which goes beyond the RKKY interaction. 

Keywords : RKKY interaction; Magnetic impurities; Spin polarized electron gas. 


1. Introduction 

The Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction between two magnetic 

moments embedded in a solid was originally derived for the free electron gas, 033 

and it has been extensively studied since then. It has the Heisenberg form JSi ■ S 2 

and shows an oscillatory behavior as a function of the distance R between the two 

magnetic moments; for example, J(R) oc (xcosx — sinx)/cc 4 for the electron gas 

in three dimensions (3D), where x = 2 kpR with kp being the Fermi momentum. 

Using lattice models, Dzyaloshinsky and Moriya (DM) showed thatEEl in situations 

with broken symmetry (inversion or time reversal), there appears a vector or a 

—> ^^ 

tensor interaction term, D and r, respectively, in addition to the scalar, RKKY 
interaction J, so that the net interaction energy is given by the Hamiltonian TL = 
J Si ■ S 2 + D ■ Si x S 2 + Si- F -S 2 . There is currently a considerable interest on this 
less-studied DM interactions, because of the recent observational] of skyrmions, 
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a novel vortex-like spin structure^, that forms due to the competition between the 
RKKY and the DM terms. 

While the RKKY interaction is quite simply illustrated by taking the example of 
the free electron gas, this is not the case for the DM interactions. The original work 
of DM involved lattice models, while in many subsequent works, complicated inter¬ 
actions such as non-magnetic spin-flip centers 10 or relativistic scattering potentialJ 11 ^ 
were used to obtain the DM interactions. In this paper, we point out that the simple 
generalization of the electron gas to include spin polarization does indeed produce 
the DM interaction, albeit just the tensor part as the vector part D = 0 because the 
inversion symmetry is still present. 



Fig. 1. (Color online) The interaction between two magnetic moments mediated by the host 
electrons in the crystal, which leads to the well known RKKY interaction of the type JS 1 ■ £ 2 - 
However, if the host crystal has broken symmetry, then additional interaction terms appear, known 
as the Dzyaloshinsky-Moriya interaction. 


Before we discuss the case of the spin-polarized electron gas, we note that the 
coupling between the spin and orbital degrees of freedom is essential to the DM 
interactions, as without it electron energies are spin degenerate, e{k f) = s(k |), 
and it turns out that, if this is true, then there would only be the RKKY interaction. 
This is the case for the standard, non-spin-polarized electron gas. If there is sufficient 
symmetry, even in the presence of spin-orbit coupling, the electron states may still 
be spin degenerate, so that the presence of the spin-orbit coupling alone is not 
enough to produce the DM interaction. 

In general, the presence of the inversion symmetry in a solid leads to the condi¬ 
tion 


£ k J [ ~ £ -fcf 

while the time reversal symmetry implies that 


( 1 ) 


so that if both symmetries are present, one obtains the spin-degenerate energies, 
viz., = £jr^. This equality is violated if there is at least one broken symmetry. 
The simplest way to include such a broken symmetry is to consider spin polarized 
bands in a 3D electron gas, which has broken time reversal symmetry. This may 
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for example be produced by a magnetic field, either externally applied or produced 
by the constituent atoms in the solid. We now consider the interaction between 
magnetic moments placed in the spin-polarized electron gas. 


2. Spin-polarized electron gas and magnetic interaction 


The magnetic interaction energy between the two localized moments is obtained 
by using the quantum mechanical perturbation theory. The first magnetic moment 
perturbs the electron gas, which is then felt by the second moment, thereby coupling 
the two moments. The unperturbed eigenstates of the electron gas are plane wave 
states 


| ka) = —!=e*k'’ 

1 ' Vn 


I o'). 


n 2 k 2 

2 m 


-F A, 


( 3 ) 


where we have taken the energies to be spin dependent with the sign correspond¬ 
ing to the energy of the spin f and states, respectively, and the box normalization 
has been used with Q being the volume of the box. We now introduce two localized 
spins S i and S 2 , located at the origin and at the position R , respectively, which 
interact with the electron gas via the contact interaction 


V 1 (r) = -XS(r)S 1 -s, 

V 2 (r) = -X5(r-R)S 2 -s, (4) 

where s is the spin of the electrons, and, in operator form, the perturbing potentials 
may be written as V\ = — A|r)(r|S'i • s and similarly for V 2 . 


Energy 



k 


Fig. 2. (Color online) Electron gas with spin dependent energy. Ep denotes the Fermi energy 
and the Fermi momenta in the two spin channels are indicated by kp^ and kp±. 


The perturbed electron states \ka)) due to the first moment may be expressed 
in terms of the Lippmann-Schwinger equation^, correct to the first order in per- 
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turbation theory 

| ka)) = \ka) + GVi\kcr), (5) 

where G(E) = (E + i[i — H ) _1 with fj, —> 0 + is the unperturbed retarded Green’s 
function to be evaluated at the energy E = £r , and H is the unperturbed Hamil¬ 
tonian that produces the eigenstates in Eq. (|3|. Eq. ([5| may be written in terms of 
the matrix elements of the Green’s function, which reads 

|fccr» = | ka) - ~^= J d 3 r \r(j')G ac7 ' (r, 0, ■ s a<7 ,. (6) 

The Green’s function in the present case is diagonal in the spin indices 


G aa '(r,f",E) = ( ra\G(E)\r'a') 

_ 1 e ik-{r-r ') 

-nZgE + i»- Hta 


( 7 ) 


The spin S 2 at R feels the perturbation via the contact interaction V 2 (r), leading 
to the interaction energy 

occ 

E(R) = ^«fca|E 2 | ka)) = E 0 + E(R) + 0(A 3 ). (8) 

ka 

This expression is evaluated by substituting Eq. |o]) for the perturbed states \ka)) 
and performing the summation over the occupied states. The Eq term, which is 
0(A), introduces a constant shift in energy; In fact there is an additional 0(A) 
term, viz., (ka\V\\ka) , originating from the perturbation of the first spin to the 
electron gas. Adding this to Eq. (|8| to redefine Eq, and performing a straightforward 
momentum summation, it is easy to show that, for the 3D electron gas, 

Eo = ^(S lz + S 2z )(k 3 Ft - k 3 n ). ( 9 ) 

A similar expression for the 2D electron gas is given later. 

This result is easily interpreted, viz., that Eq is the product of the energy gained 
per electron £0 times the excess number N of spin up electrons over spin down 
electrons, where £0 = —(Afi/211) x (Sis + S 2z ), as computed for a single plane- 
wave electron using the first-order perturbation theory with the perturbation Eq. 
Q, and N = (47r/3 ){k F ^ — kp^)/(8ir 3 /fl). Note that we have used the standard 

Pauli matrices for the itinerant electron spins, e.g., s z = ( H/2 ) Q , etc., so 

that effectively we have chosen z as the spin quantization axis and this is why the 
z components S\ z and S 2z appear in Eq. ©• 

The next term in Eq. (|8|, which is 0(A 2 ), produces the magnetic interaction 
and may be written as 

x 2 occ 

E(R) = ^ ^(a \S 2 ■ sG(R, 0,£ fc ) S, ■ s + c.c. (10) 

ka 
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Now, as mentioned already, G G ' a " ( R , 0, ) is diagonal in the spin space, though it 
can be different for a =f or for a given k. It is convenient to write them in terms 
of the Pauli matrices: 


G(R, 0, £g t ) = g 0 cr 0 + -u{<jq 


G(R, 0, £g,) = 3ocr 0 + -d(a 0 + a z ), 


(ID 


-kiJ • 2' 

where go, u, and d are c-numbers and functions of k and R. Putting this into the 


energy expression Eq. (101, we now get 
A 2 ( occ 

E(R) = VV ^(<r|S 2 ■ s g 0 a Q Si • 8 \a)e~ tk ’ R 


ka 


+ l^ 2 ' * ^ 9 ° a ° + |(o- 0 - a z ))Si ■ s | f)e lk R 


k =0 
hF\, 


+ ^4 |S 2 • s(g 0 a 0 + -(cr 0 + cr 2 ))Si • s | 4_)e lk R + c.c. 
k -o ) 

This can be evaluated using the spin identities 
Tr[Si ■ a S 2 ■ a] = 2Si • S 2 , 

(t |Si • a S 2 ■ 3 | t) = Si ■ S 2 +*(& x S 2 ) z , 

4 |Si • 5 S 2 ■ 3 | = St ■ S 2 - i(& x S 2 ) z , 

(t |Si • 3 a z S 2 ■ 5 11) = ~Sx ■ S 2 - i(S i x S 2 ) z + 2 S lz S 2z , 

4 |S! • 8 a z S 2 ■ 5 U) = S x • S 2 - i[S y x S 2 ) z - 2S 1z S 2 ,, 


and the energy expression Eq. (12) can be cast into the desired form 


E(R) = JS 1 - S 2 + D ■ (Si x s 2 ) + Si • r • s 2 , 


( 12 ) 


(13) 


(14) 


where the various coefficients may be obtained by converting the summations into 


integrations in Eq. (12). 


3. Results 

3.1. Electron gas in 3D 

We now evaluate the interaction energy coefficients in Eq. ( |14[ ) for the spin polarized 
electron gas in 3D by evaluating the Green’s function and performing the momentum 
integrations. A straightforward contour integration 13 4f the expression ([T]) yields the 
Green’s function 

0 ia(E+A)r 


G(r, 0, E) = 


—2 m ( e 


AnrH 2 


0 


0 

p ia(E—A)r J ’ 


(15) 
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where a(x) = (2 mh 2 x ) 1 ^ 2 for x > 0 and i(2mh 2 |cc|) 1//2 for x < 0. Comparing Eq. 

@_ 

> A, and d = p e ^ k ' 2+2S ' 2R - 


15|) with Eq. (|11|), we have: go = /3e lkR , u = (3e v /2<52 k2R — g 0; if _A < < A, 

> A. Here 


0 Wk 2 -25 2 R _ 


if £ r 


or /3e Wk 2d H - ffo , if £g f > a, and a = pe"*" 1 "" “ - g 0 , - 
P = — m(2Trh 2 r)~ 1 . Performing the necessary integrations, we finally get from Eqs. 
fl2|) and (fl 


the result for the RKKY interaction 


J = - 


A 2 m 
(47r) 3 i? 4 


x ( J Ot) + ^(a:;)), 


(16) 


where x-j- = kp^R, = kp^R, and I(x) = sin(2x) — 2xcos(2x). Interestingly, the 
vector DM interaction D = 0, while the tensor interaction survives: 


(P 0 0\ 

r= o p o , (17) 

\o o 0/ 


where P = g{x±, X|) — J(x±, £f) and defining g = (x 2 — x 2 ) 1//2 = y/25R, the expres¬ 
sion for g is 


g(xi,x t) = - 


m 


8t r 3 R 4 

XI 


cos 


lo 


o-V 7 ! 2 -* 


r x t 


cos 


\Jx 2 — g 2 


\/x 2 + g 2 ) x xsinx dx. 


(18) 


We note that for the unpolarized electron gas, putting kFf = A'fj, = kp, the 
above expressions yield P = 0, so that we recover the standard RKKY result: 
E(R) = J Si ■ S 2 , where J = —\ 2 m./(32ir 3 R 4 )I(kpR). We illustrate the oscillatory 
behavior of J and P in Fig. ([3]) for typical parameters. The beating pattern seen 
in the figure comes from the interference of the two Fermi momenta for the two 


spins. For example, it is easy to see from Eq. (161 that for large distances J ~ 
R~ 3 cos(kF^R + kp^R) x cos(kp^R — kp^R), which leads to the beating pattern as 
indicated in Fig. Q. 

Another point is that the constant energy shift Eq, Eq. |9]), is proportional to 
XSkF , where Skp = kp^ — kp^ is the spin polarization at the Fermi surface, while the 


magnetic interactions, Eqs. (16 


18), are proportional to the second order A 2 . Thus, 
the latter terms dominate for the weakly spin-polarized case, viz., Skp « (mkp)A, 
where kp = ( kp f + kp^)/2 is the spin-averaged Fermi momentum and the typical 
scenario, Skp « kp, has been assumed. Finally, note from the form of the tensor 
interaction Eq. ( p~7] > that the net interaction Eq. (14) can be written as a sum of the 
Heisenberg and Ising interactions, E(R) = J' S\ ■ S 2 + J"Si z S 2z , where J' = J + P 
and J" = —P. This is true both for the 3D electron gas as well as for the 2D gas 
described below. 
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Fig. 3. (Color online) The RKKY interaction J (top) and the tensor Dzyaloshinsky-Moriya in¬ 
teraction P (bottom), in units of mn~ 3 X 2 , for the 3D electron gas as a function of distance R, 
obtained from Eqs. ( |l6| - |18) . 


3.2. Electron Gas in 2D 

In 2D, the Green’s function Eq. 0 becomes 


G a o'{r,r',E) = 


p ik-(r—r ) 


(27r) 2 J E + i/i — et 


d 2 k. 


(19) 


Using the Jacobi-Anger expansion 1 ^ in terms of the Bessel’s functions, viz., e lk ' r = 

OO 

Jo(kr) + 2 J2 i n Jn(kr) cos (nQ), where 9 is measured with respect to f, the angu- 

n—1 

lar integration produces just one non-zero term involving Jo(kr), and the radial 
integration finally leads to the well-known resultP^ 


G n 


r (P,f',E) = -^K 0 



f'|(A±A) 1/2 , 


( 20 ) 


where + (-) is for er =f (J,) and K 0 is the modified Bessel function of the second 
kind. For the perturbation of the ka state, as before, we put E = er and express 


Eq. (201 as a function of the momentum k. The result is: 


G(R, 0, 


■kt 


\ _ or (K 0 (-ikR) 0 

> P \ 0 K 0 {-i{k 2 -2S 2 ) 1 / 2 R) 


( 21 ) 


where 6 = (2mA) 1 ^ 2 H 1 and f}' = —m/(TrH 2 ). The same result for the down spin is: 


G(RAs h ) 


, (K 0 (-i{k 2 +25 2 ) 1 / 2 R) 0 \ 

P ( 0 Ko(-ikR))' 


( 22 ) 


Writing this in terms of the Pauli matrices as in Eq. (JTTj) , we get: go = 
(3 1 K 0 (—ikR), u = /3'K 0 (V^S 2 — k 2 R) — go, if —A < < A, and (3 1 K 0 

(— iy/k 2 — 2 S 2 R) — go, if > A, and d = /3 'Kq (— iy/k 2 + 2 S 2 R) — go, if > A. 
The remaining algebra is the same as the 3D case and the final results for the 2D 
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electron gas reads: 


E o° = ( S u + S 2z){xf - x\), 

j2D= X ( / 2( a: l) + / 2( x t)) 5 

D 2D = 0 , 


(23) 

(24) 

(25) 


p2D _ 


m\ 2 

8nR 2 


x - 


[ Re[K 0 (y/rf - x 2 )] + [ Y 0 (^x 2 - rf) 

J 0 j 71 


+ J Y 0 (\/x 2 + r? 2 )^ x xJq(x) dx - J 2D . 


(26) 


where, again, rf = 2 S 2 R 2 1 and / 2 (a:) = x 2 [J 0 (x)Y 0 (x) + Ji{x)Y 1 {x)\, the J’s and F’s 
here being the Bessel’s and Neumann’s functions, respectively. In the large distance 
limit, ^(x) —> —7T _1 sin(2a;), so that the magnetic interactions fall off as ~ R~ 2 , 
with a beating pattern similar to the 3D case, as might be seen from Eq. (251 
for J 2D . Once again, in the non-spin-polarized limit, kp-[ = kp ±, we recover the 
standard results for the RKKY interaction for the 2D electron gas. 


4. Summary 

In summary, we extended the original RKKY description for the magnetic interac¬ 
tion in the electron gas to the spin-polarized case. The broken time-reversal sym¬ 
metry produces, in addition to the RKKY interactions, the DM interaction terms, 
explicit expressions for which were derived for the cases of 2D and the 3D electron 
gas. Numerical results showed a beating pattern due to the interference between 
the two Fermi momenta kp ^ and for the two different spin channels. In the 
non-spin-polarized limit, kp^ = kp±, the DM interaction vanishes and the standard 
result for the RKKY interaction J for the electron gas is recovered. Since for the 
spin-polarized electron gas, treated here, the inversion symmetry is not broken, the 
vector DM interaction term vanishes. We note that in opposite situations where the 
inversion symmetry is broken, e.g., for the electron gas with Rashba or Dresselhaus 
terms, the vector DM interaction D is non-zero. 1 ZiIH 
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